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In a recent work [Phys. Rev. Lett. 108 , 136401 (2012)] we have addressed the problem of a 
Luttinger liquid with a scatterer that allows for both coherent and incoherent scattering channels. 
We have found that the physics associated with this model is qualitatively different from the elastic 
impurity setup analyzed by Kane and Fisher, and from the inelastic scattering scenario studied 
by Furusaki and Matveev, thus proposing a new paradigmatic picture of Luttinger liquid with an 
impurity. Here we present an extensive study of the renormalization group flows for this problem, 
the fixed point landscape, and scaling near those fixed points. Our analysis is non-perturbative in 
the elastic tunneling amplitudes, employing an instanton calculation in one or two of the available 
elastic tunneling channels. Our analysis accounts for non-trivial Klein factors, which represent 
anyonic or fermionic statistics. These Klein factors need to be taken into account due to the fact 
that higher order tunneling processes take place. In particular we find a stable fixed point, where an 
incoming current is split \ between a forward and a backward scattered beams. This intermediate 
fixed point, between complete backscattering and full forward scattering, is stable for the Luttinger 
parameter g < 1. 


I. INTRODUCTION 

The concept of a Luttinger liquid (LL) provides a very 
general framework to deal with a strongly interacting 
electron gas confined to one spatial dimension (ID) 1 . 
Contrary to higher dimensional situations, where the 
quasi-particle (qp) concept describes low energy excita¬ 
tions, in a LL only collective excitations are long lived 
and stable. Particle-like excitations of a LL have an en¬ 
ergy dependent density of states and can relax in en¬ 
ergy in the presence of backscattering. Applications of 
the LL concept include semiconductor 7 , metallicP, and 
polymei® nano wires, carbon nanotube^, quantum Hall 
edgetP, and cold atom^. 

Due to the geometrical restriction to ID, electrons may 
be divided into two sectors, right moving and left mov¬ 
ing. The presence of an impurity®!gives rise to inter¬ 
branch scattering, i.e. backscattering with a power-law 
dependence of the scattering probability on energy. Hith¬ 
erto there have been two paradigmatic models which ad¬ 
dressed impurity scattering in the context of LL: a purely 
elastic impurity as discussed by Kane and Fishei® (KF), 
and totally inelastic scattering described by Furusaki and 
Matvee^ (FM). 

In a recent publication 11 -! we have introduced and stud¬ 
ied a model which is a hybrid between the two. We 
have addressed interacting electrons in one dimension de¬ 
scribed by a Luttinger model, which includes a scatterer 
that may give rise to elastic coherent tunneling, and at 
the same time may accommodate inelastic modes. Our 
analysis has led to predictions which are qualitatively 
different from the KF and the FM pictures, and in this 
sense can be considered as a new paradigmatic scheme of 
impurity scattering in the context of LL. 

In particular, we have obtained the following results: i) 


Due to the presence of both elastic and inelastic scatter¬ 
ing channels, there exists a new stable non Fermi liquid 
fixed point (FP). Asymptotically, the impurity becomes 
a symmetric beam splitter, which halves the incoming 
beam into two equal outgoing beams. For this reason 
we termed this limit a | - | FP. We also identify other 
non Fermi liquid FPs, which are stable in a certain di¬ 
rection and unstable in another direction in parameter 
space. A major facet of our earlier work was that upon 
renormalizing our model down in temperature (or bias 
voltage), under generic conditions the model flows to a 
stable FP. The neighborhood of this FP is marked by 
the non-Fermi liquid correlations described above, ii) 
At equilibrium, there is thermal noise at each incoming 
and each outgoing terminal, but no cross-correlation in 
the thermal noise, i.e. no correlations between incoming¬ 
incoming and outgoing-outgoing edges - in stark contrast 
to the Landauer-Buttikcr picture, iii) Out of equilibrium, 
when the system is voltage biased at one of the incom¬ 
ing terminals, the impurity acts as a \ - ^ beam splitter 
with no shot noise component in the outgoing current. 
Similarly, cross-current correlators do not contain a shot 
noise component. 

Our analysis in reference m involved bosonization of 
the interacting fermionic Hamiltonian. It was based on 
perturbation in the elastic (coherent) scattering ampli¬ 
tudes in the presence of a charging term representing 
electrostatic correlations in the scattering region. We 
have concluded that for energies below the charging en¬ 
ergy, the bosonic fields affected by the elastic scattering 
terms become massive. As a result, only one of the four 
independent bosonic fields in the model remains asymp¬ 
totically free. The asysmptotically Gaussian action has 
facilitated calculation of noise and current-current corre¬ 
lations. 
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Naturally, the fact that our analysis has been based 
on, and the results have been motivated by perturbative 
analysis, raises some questions. Most important is the 
fact that various relevant terms of the action have been 
analyzed within a renormalization group (RG) scheme 
separately. Evidently this procedure ceases to be justi¬ 
fied when the respective amplitudes of the elastic tun¬ 
neling terms grow through the RG analysis. One may 
need to evaluate expectation values of products of non¬ 
commuting terms, which a naive perturbative analysis is 
incapable of doing. Careful analysis is required in this 
process. 

We note that the elastic tunneling channels involve 
scattering from any of the two incoming channels to any 
of the outgoing channels. To establish the right language, 
and to enable a compact description of the RG flow di¬ 
agram, we divide the elastic scattering terms into back 
scattering (the terms connecting channels 1 with 3, and 
2 with 4, cf. Fig. 2) and forward scattering (the terms 
connecting 1 with 4, and 2 with 3). We assume that the 
coefficients of the two back scattering terms are equal 
to each other, and similarly the coefficients of the two 
forward scattering are equal to each other. Relaxing this 
assumption does not modify the asymptotics of the prob¬ 
lem, but makes the description of the RG flows more cum¬ 
bersome. Assuming we have renormalized the problem 
down from high temperature (or voltage), and that cer¬ 
tain irrelevant non-universal terms have been renormal¬ 
ized out, we are left with flows in a two-dimensional pa¬ 
rameter space. We indeed establish the existence of a new 
stable \ - \ FP, which corresponds to strong coupling 
in both scattering channels (backscattering and forward¬ 
scattering respectively). We denote this a strong-strong 
fixed point (SSFP) . In addition, we find and character¬ 
ize three more FPs: a trivial unstable weak scattering FP 
(weak-weak fixed point, WWFP), and two other FPs: a 
weak-strong fixed point (WSFP), and a strong-weak fixed 
point (SWFP), each being stable in one direction of pa¬ 
rameter space and unstable in the other direction. All 
these points are marked by non Fermi liquid correlations. 

We note that the issue of non-commutativity of the 
various tunneling operators is central to the analysis of 
the model beyond the perturbative limit. Here, we will 
be dealing with operators describing tunneling of Abelian 
quasi-particles. Such particles (known as anyons) possess 
fractional statistics, intermediate between bosons and 
fermions. The exchange of two qps involves a statisti¬ 
cal phase which is ill-defined: its sign depends on the 
’’history” of the trajectories employed in the course of 
the exchange. In order to avoid such a pathology, one 
may resort to either one of the three following tricks: i) 
avoid single qp field operators (qp creation or annihila¬ 
tion operators) and assign statistical Klein factors only 
to operators which are bil inea r in qp field operators, e.g. 
to qp tunneling operator :) 13 f |4 ( ii) attach a statistical flux 
tube to a qp an d follo w the kinetics by way of a quantum 
Master equatiorl— iii) introduce all relevant edges on a 
single contour, and choose a convention how qp’s are be¬ 


ing exchanged 1 -^ (i.e., clockwise or anticlockwise). In the 
present analysis we will assign generalized Klein factors 
to qp operators, following the philosophy of iii). 

The outline of the paper is the following: In Section 
II, we introduce the problem addressed here, and discuss 
its mapping onto a model setup which allows system¬ 
atic bosonization and subsequent analysis. The assump¬ 
tions underlining this modelization and the applicability 
thereof are outlined. In Section III, we position our build¬ 
ing blocks (in particular the chiral channels of our model) 
within a particular geometry, giving rise to a well-defined 
convention that determines the commutation relations 
of the anyonic quesi-particles. Performing a sequence of 
canonical transformations allows us to derive an effec¬ 
tive action that captures the symmetries of the problem. 
Section IV is devoted to the perturbative analysis of the 
weak scattering fixed point (WWFP), while Section V is 
focused on the study of the scaling near the weak-strong 
fixed point WSFP. We note that the physics of the strong- 
weak fixed point SWFP is trivially obtained by trivial 
exchange of field indices from the WSFP. Here, we argue 
that the low energy dynamics is dominated by a phase 
slip-instanton mechanism. In Section VI, we analyze the 
scaling near the SSFP. In this limit, an instanton picture 
within a 2-dimenional parameter space is employed. In 
the limit of non-interacting Luttinger wires our analysis 
needs special care: some of the flows become marginal. 
This is discussed in Section VII, where comparison with 
known results is presented as a consistency check. Fi¬ 
nally, in Section VIII, we present a short summary of our 
results, and a few proposals and speculations concerning 
further directions of our analysis. In Appendix A, we 
make a few comments concerning experimental verifica¬ 
tion of our predictions, and in Appendix B, we compare 
the results obtained within our bosonic theory with an 
exact refermionization treatment for a specific value of 
the Luttinger parameter. 


II. THE SETUP AND ITS MODELING 

Our model is motivated by an experimental setup with 
a two-dimensional electron gas in a strong magnetic field, 
such that the fractional quantum Hall (FQH) regime is 
reached. We focus on a Laughlin filling fraction, such 
that, in the absence of edge reconstruction, the structure 
of the edge is simple, consisting of a single chiral channel 
at each edge. In addition, we imagine the existence of 
a compressible puddle with somewhat reduced density, 
e.g. v = 1/3 in the bulk and a puddle with a compressible 
v = 1/4 state of composite fermions. Such a setup can 
be realized by locally modifying the backgate voltage in 
a certain region, see Fig. [T] 

Due to its small geometric size, the puddle has a charg¬ 
ing energy associated with it and constitutes a quantum 
dot (QD). The edge modes support anyonic qps, which 
may be backscattered through the incompressible bulk 
near the entrance and exit of the QD. The setup we thus 
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FIG. 1. Sketch of a QH bar with a gate defined QPC. The 
bulk of the Hall bar is incompressible with filling fraction 
v, the density in the QPC region is lower, such that mixing 
between edge states in a compressible region can occur, both 
via elastic and inelastic channels. 



FIG. 2. (Color online) Schematic model of a quantum 
wire with two incoming (labeled 1,2) and two outgoing chi¬ 
ral modes (labeled 3,4) connecting to an extended scattering 
region. The dark shading represents a region of capacitive 
charging, the quantum dot. Near the entrance and exit to 
the quantum dot, where translational invariance is broken, 
backscattering of quasi-particles across the Hall bar is possi¬ 
ble, denoted by vertical dashed lines. In a possible realization 
with a compressible quantum dot region, the edge current 
spreads into the bulk, and both elastic transmission along 
the edge (indicated by horizontal dashed line) and inelastic 
transmission through the dot is possible. 

envision is that of a QD with four semi-infinite chiral 
edges coupled to it. We assume that the contacts be¬ 
tween the leads and the QD are ideal, i.e. reflectionless. 
In total, one needs to consider a Hamiltonian which de¬ 
scribes the ideal chiral edges, the contact between edges 
and puddle, the puddle modeled as a metallic quantum 
dot (QD) with charging energy, and elastic (coherent) 
tunneling channels that may support qp tunneling (see 
Fig. 2). In terms of modeling this setup, we need to find 
a way to describe an ideal, reflectionless contact between 
a chiral lead and a QD. In order to achieve this goal, 
we follow Matveev: we employ four electrostatically cou¬ 
pled infinite chirals. For each of these infinite chirals, 
one half is identified with the chiral wire outside the dot, 
whereas the other half is assumed to be positioned inside 
the dot. In this way, low energy excitations of the chiral 
parts ’’inside the dot” can be identified with dissipative 
processes within the QD. Transport processes between a 
chiral entering the dot and another chiral leaving the dot 
are mediated by i) the charging energy within the dot, 
which prevents charge accumulation inside the QD, and 
ii) by coherent elastic tunneling between incoming and 
outgoing chirals. The setup corresponding to our model 
is depicted in Fig. 3. 

We are now in a position to write the total Hamiltonian 


. It consists of three terms: The low energy dynamics of 
the original semi-infinite leads, as well as the part of the 
chiral edge channels that constitutes partof the quantum 
dot, are delegated ot a Hamiltonian that represents four 
infinite chirals, iT c h- Separately, the charging energy of 
the QD is described by H dot . Finally, the coherent part 
of the lead-QD tunneling is represented by H tun . We re¬ 
sort to bosonization, and define bosonic fields (f>i ( x ) with 
i = 1,2,3,4 for each chrial, see Fig. 2. We decompose 
the bosonic fields into finite momentum and zero modes 
according to 

<fii{x) = <j>i(x) + $(x) . (1) 


The finite momentum part has a Fourier representation 
according to 



and the structure of the zero modes will be discussed in 
the context of Klein factors in Section III. For the mode 
decomposition above, the plus sign should be used for 
the right-moving modes 1 and 4, and the minus sign for 
the left moving modes 2 and 3. The and a\ are 
canonical boson operators with commutation relations 
[ a i,q’ a j,q'} = 27T<5(qr—qr'), giving rise to [<Pi{x),(pj(y)\ = 

(Trg/2)sgn(x — y). The charge density is related to be bo¬ 
son fields via pi(x) = x), and the Hamiltonian 

for the chirals is given by 

H ( ;h = Y, f dx [d x <pi(x)] 2 . (3) 

11 y i—i •* 

In addition, there is a Coulomb charging energy for the 
dot 

H dot = ^Q 2 , (4) 

where 

Q = - 7 =- [-01 (0) - 0 2 (O) + 03(0) + 04(0)] (5) 

denotes the charge on the dot. Coherent elastic tunneling 
between incoming and outgoing wires is described by the 
Hamiltonian 

-fftun = t\3 COS (if! - lfi 3 ) + t 2 4 COS (ip 2 - <£4) ( 6 ) 

-K14 COS (ipi — cp d ) + t 2 3 COS {tp 2 — p 3 ) ■ 

Here we ignore retardation effects, assuming that tunel- 
ing is instantaneous. Thus, the total Hamiltonian for the 
system is given by 

H = H ch + H dot + H tun . (7) 

As for some parts of our analysis it will be useful to resort 
to a functional integral formulation of the problem, we 
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FIG. 3. (Color online) Effective modelling of the system (cf. 
Fig.[2f through four infinite chiral channels. Each channel 
consists of a half infinite wire (either incoming or outgoing), 
and a half infinite wire which is part of the quantum dot 
(shaded area). Elastic scattering channels are represented by 
wiggly lines. 


also state the Matsubara finite temperature action of the 
system 



Here, the plus sign refers to the right moving branches 1 
and 4, and the minus sign to the left moving branches 2 
and 3. We note that the zero mode part of the Hamil¬ 
tonian H 0 and and of the action Sq will be discussed in 
detail in Section III. 

A few commons are in order now: (i) Our analysis will 
be performed using the Matsubara technique, (ii) The 
tunneling is that of anyons (qps). We consider only this 
type of tunneling terms as they are (when allowed) po¬ 
tentially the most relevant terms, (iii) All channels are 
put on an equal footing. There is no distinction between 
forward scattering and backscattering any more. Par¬ 
ticles on different chirals are different species, hence the 
finite momentum boson fields <f> i, <f> 2 , commute. It is 
therefore vital (for the sake of preserving the appropriate 
commutation relations) to introduce Klein factors which, 
in the case of qp tunneling, will be anyonic Klein factors. 
In the notation of Eq. the Klein factors are implicit 
in the prefactors tij, their precise realization will be dis¬ 
cussed below. We denote the c-number valued scattering 
amplitudes without Klein factors by 7 ^, i.e. tij = 'Jij'x 
(Klein factors), (iv) Each infinite chiral contributes half 
to the lead (either incoming or outgoing) and half to the 
degrees of freedom of a compressible puddle. The puddle 
(a.k.a. quantum dot) accommodates gapless soft modes, 
(v) The charging (that is, electrostatic) interaction intro¬ 
duces coupling among the 4 channels. As we renormal¬ 
ize the model down in temperature (or applied voltage), 
flowing towards low energy dynamics (long time scales), 
the presence of a charging energy scale will force charge 
fluctuations on the QD to freeze out. The QD then re¬ 
mains effectively neutral, (vi) The model is a hybrid of 
elastic (coherent) scattering between incoming and out- 
goingleads, and inelastic transport through the QD. The 
latter involves the excitation of low-lying modes in the 


QD, which can be described in terms of many particle- 
hole excitations. The former represents (to lowest or¬ 
der in the tunneling coupling constant) elastic tunneling 
through the QD. While one may imagine that backward 
tunneling, say between channels 1 and 3, is mediated by 
a resonant impurity level in a Hall bar constriction, one 
can motivate the forward tunneling, say between chan¬ 
nels 1 and 4, as being due to coherent propagation along 
the edge of the QD. For any real QD, there is an energy 
scale defined by the mean single-particle level spacing in 
the QD, A. Once the maximum of temperature and bias 
voltage is driven below A, the soft modes which consti¬ 
tute an essential part of our model are frozen out. At that 
point our model ceases to be valid. Nevertheless, as long 
as the infrared cutoff of our RG flows is larger or compa¬ 
rable to A (which might be extremely small) the model 
and the predictions following from its analysis, are valid, 
(vii) The interesting regime of our RG flows stipulates 
running frequencies A <C |w| E c . Modes involving 
charge fluctuations are frozen out; for the model to be 
valid, soft modes in the QD still need to be available. Is 
such a regime possible? If we are focusing on a QD sub¬ 
ject to a strong perpendicular magnetic field, then the 
relevant spectrum of the QD is Id (on the edge). Osten¬ 
sibly, both E c and A will scale like 1/L (L= the linear 
size of the QD), and the difference between them may 
be only due to different coefficients. However, we may 
control the charging energy E c , and actually render it 
larger than A, by considering a puddle realized by a sea 
of composite fermions, for instance at filling v = 1/4, 
whose motion does not follow Id edges but instead ex¬ 
plores the available 2d area of the QD. 

Our analysis presented here is done with a view to the 
chiral edges of a FQHE. We stress here that it may apply 
to quantum wires described by Luttinger liquid models. 
Applying our analysis to Luttinger liquids implies certain 
assumptions. This is discussed in Appendix A. 


III. GENERALIZED KLEIN FACTORS, AND 
THE EFFECTIVE ACTION 


As was commented in Section II, we need to establish a 
convention as far as Klein factors, underlying the anyonic 
field operators, are concerned. We choose to establish 
such generalized Klein factors in relation to a geometry 
where all 4 infinite chirals are assumed to be segments of 
a single chiral contour. In that convention, the exchange 
of two anyons positioned on channel i and channel j re¬ 
spectively has an unambigious chirality. The geometry is 
shown in Fig. 4, and the related bosnonic field configu¬ 
ration is depicted in Fig. 5. 


A. Effective Geometry 

As we will consider qp scattering between the wires, we 
need to find a faithful representation of quasiparticle op- 
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FIG. 4. Putting all four chiral wires onto a single contour. 
Our previously disconnected chiral wires become segments de¬ 
noted 1,2,3,4. They are connected by scattering as shown (cf. 
Eq. Q. The quantum dot is represented by the shaded area. 
All segments should be imagined ’infinitely long’, e.g. a quasi¬ 
particle entering the dot along 1 will then propagate vertically 
upwards but will never make it to the ’top of the contour’ in 
finite time. It has to scatter to get to 4 (henceforth to 2). Note 
that the vertical scattering amplitude 7 „ = I 713 + e _I7r ® 724 
and the horizontal scattering amplitude 'yn = I 723 +7m| will 
be the parameters of our RG flow as shown in Fig. [0 The rel¬ 
ative phase factor in the definition of y v is due to the phases 
of scattering operators in Eq. (171. 


erators in terms of bosonic constituents. For the purpose 
of establishing unique commutation relations between qp 
operators acting on different chirals, we will put the four 
chiral edges defining our problem onto a single contour, 
as shown in Fig. [4j We now follow the recipe for repre¬ 
senting quasiparticle operators as discussed by Oshikawa, 
Chamon, Afleck (OCA in the following) in Appendix E5 
of Ref. [IT]. The starting point is a representation of 
chiral bosonic field, generalized for the presence of zero 
modes. In operator notation, 


<Pi{Xi) = -J= ; PiXiS , 




(9) 


Here, 27 is confined to the jth segment of the contour, and 
L is extension of that contour. The sign factor is S{ = 1 
for i = 1,2 and s,; = —1 for i = 3,4. The Fourier compo¬ 
nents of the oscillator modes are such that these modes 
are algebraically independent on the different contours, 
and hence they commute. The role of zero modes and 



FIG. 5. Field configuration on the closed contour. We intro¬ 
duce zero modes qi, pi to describe the increment of the field 
throughout any of the four segments. The ’finite momentum 
components’ 4>i obey periodic boundary conditions. Shaded 
regions indicate the dot regions. 



FIG. 6 . Description of the RG flow in terms of masses 7 h and 
7 „, which denote the strength of elating tunnelling in horizon¬ 
tal and vertical direction, respectively. Near the RG-unstable 
weak-weak fixed point 1 , there are perturbative corrections 
to the Gaussian (1/2,1/2) conductance in both yu and y v , 
which grow under the RG. In the vicinity of the weak-strong 
fixed point 2 , instantons in the strong scattering term y v are 
irrelevant, constituting a stable direction in parameter space, 
while the weak scattering amplitude yn grows under the RG 
and thus constitutes an unstable direction. The conductance 
is close to unity in the vertical direction, and perturbatively 
small in the horizontal direction. In the vicinity of the strong- 
strong fixed point 3, both vertical and horizontal tunnelings 
are strong, with RG irrelevant instanton corrections. The 
conductance is ( 1 / 2 , 1 / 2 ), and the fixed point is stable under 
RG scaling. 


finite momentum modes in manufacturing a field contin¬ 
uously covering the whole contour is indicated in Fig. [5] 
We note that the zero modes satisfy the commutator re¬ 
lation 


[Qj,Pk] = giSj,k- ( 10 ) 

At this point, a comment on normalization is due: we rep¬ 
resent our chiral fields as exp(iv / 2 v 5 ) ~> exp[i (<7 — yyP x + 
y/2<j>)]. This differs from the OCA convention by the fac- 
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tor \/l/g. OCA have a convention exp [i(y/gq+ 27T ^ p + 
(j>)], and unit commutation relations [q,p] = i. 

When substituted into a functional integral, qi qi (t) 
becomes a time dependent field. In the literature, one 
often finds the additional dynamical substitution rule, 
PiXi —> pi(xi~t), where, in order to simplify the notation, 
we have taken the velocity to be unity. This we should 
interpret in terms of the interaction picture dynamics 
generated by the zero mode Hamiltonian 


Hr 


— 71 - 2 

~ 4 gL^ Vi 


( 11 ) 


and the corresponding zero mode action 


So[q,p] 



(^ \id T q, - H 0 


The corresponding equation of motion 

dtqi = jPu 
dtPi = 0, 


( 12 ) 


define new operators T, = rjiipi, which are designed such 
that 

T, : d-v = TjVl/, _ (16) 

The derivation of these commutation relations from 
the current path integral formalism is detailed in Ap¬ 
pendix [A] 

We aim to understand how the presence of these ’Klein 
factors’ will interfere with the scaling of the scattering 
operators. To achieve this, we define the latter in such a 
way that the tunneling points sit at Xi = 0. In order to 
avoid ambiguity, we actually remove the tunneling sites 
a little bit from zero as indicated in Fig. [4] This means 
that the piXi terms in the quasiparticle amplitudes drop 
out. We are then led to consider the following bilinears: 

vf 3 vl/ 1 = g*(93-gi) e *f (j>i+P3+2p2+2p 4 ) e *(<?h-03) e -*f a ^ 

vf 4 vJ/ 2 = e i («4-g2) e *f (-P2-p4)gi(</’2-<M gif 9 

vj/ 3 Vl/ 2 = e *(93 —92) e »f (P3+P2)gi(02-03) g-if 9 

_ e *(g4-gi)g*f (pi+P4) e *(0i-04) e -if a 

(17) 


are solved by qi(t) = qi + Pij^t, with constant ( qi,Pi ). 
Substitution into the field definition Eq. © leads to the 
dynamical interpretation 




^ (qi - \pi{x iSi - t)j + <t>i{xi,t) 


In deriving these relations, we made use of the commu¬ 
tation relations 

e is f Pigis'ai _ gjs'qigisfPigiss'f g (qg) 

where s, s' = ±1. These bilinears enter the problem 
through the scattering action 


with time independent ( qi,Pi ). At this stage, the mean¬ 
ing of the sign factors s, becomes transparent. They 
make sure that on contours 1,2 causality is towards in¬ 
creasing coordinates (’incoming’), and on contours 3,4 
towards decreasing coordinates (’outgoing’). 


B. Generalized Klein factors 

Following OCA, we represent a quasiparticle operator 
obeying proper intra wire commutation relations as 


ipi = e iqi e~ i? S- piXi e^, 


(13) 


as a product of zero mode operators and the finite mo¬ 
mentum vertex operators exp(z^). Quasiparticle in¬ 
terwire commutation relations are now generated by mul¬ 
tiplication with a factor 


Vi = e l 2 ^' 
where the matrix a is given by 


oPi 


(14) 



(15) 


S a [*,*] = S v [*,*]+S h [*,*], (19) 

Sv[^, V I'] = J dr (7 13^3^1 + 7244 / 44 , 2 + h.c) , 

S h [V, T] = y dr (723^3^2 + 714^4^1 + h.c) . 

Another observable relevant to our discussion is the 
charge on the dot. Defining the charge density as 
Sid x pi = \f2rrpi, we get for the charge 
/•L 


Qi = 


1 


dx Sid x ipi = -r=—</>i(0), (20) 

2 V 27T 


Jo 

Abbreviating as </>*(0) = <f>i, we thus obtain 

Q = -\{pi +P 2 +P 0 +PA ) - ^|^(01 + 02 ~ 03 ~ 04 ) 

( 21 ) 


for the charge on the dot. (The mnemonic behind the 
different sign configuration is that in <p ~ px + <j>, the 
zero mode p comes with a spatial coordinate. That one 
changes sign when we effectively count x —> —x on the 
outgoing wires.). 

We finally note that away from the scattering points X{ 
the action is quadratic, implying that the field amplitudes 
<j>i{Xi 7^ 0) can be integrated out. As a result, we obtain 
a standard dissipative action 


The matrix structure encodes the sequential ordering of 
segments on the contour. For example, the first row in¬ 
dicates 1 < 2,3,4. The second 2 > 1,4 and 2 < 3. We 


^diss[0] — ^ ^ ' |cu n | \<fii r 

Tig 2 —' 2 —' 


( 22 ) 
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C. Effective action 


Eqs. @, and ( |B.15[ ) describe the basic con¬ 

stituents of the problem, the scattering operators, the 
charge operator, and the dissipative action in terms of the 
twelve fields and pi- Throughout, we will assume 

that the charging energy E c = e 2 /C is the largest energy 
scale in the problem. Focusing on smaller scales, we may 
assume fluctuations of the charge Q to be strongly sup¬ 
pressed, which removes one degree of freedom from the 
problem. Assuming Q to be locked, eight more linear 
combinations of the above variables may be integrated 
out without further approximation, by a procedure de¬ 
tailed in Appendix [B] To define the resulting effective 
action, we define the new set of variables, 


= 2(^1 + + fa + <M> 

$1 = 2 ~ ^2 + <^3 — 

$2 = 2 ( 0 1 ~ ^2 — 4>3 + <M> 


$3 — 2 ( 0 ! + 02 — 03 ~ <t> 4 )- 


(23) 


$3 oc Q couples to the Coulomb charging action and 
is gapped out at energy scales lower than the charging 
energy, while <3 >q is a zero mode with free action 


T 


&[<M = ^2 IWnlIdJo.ri 


(24) 


The remaining two fields <l>i j2 are the nontrivial degrees 
of freedom of the model and their effective action is given 
by 


S[$i,$ 2 ] = <Sb[i>i,$ 2 ] + S„[$2] + S^[$i], 


So[*1,*2] = 2 

y n 

S v [$ 2 ] = 7v J drcos(v / 2 $ 2 ), 

Sh[®i\ = lh J drcos(\/2$i). (25) 

Here, 7 „ = I 713 + e~™ a 7 24 | and 7 ^ = | 7 23 + 714 1 are 
effective scattering amplitudes characterizing the back¬ 
ward and forward strength, respectively, cf. Eq. (B.18). 


T (\U n \ -U r , 

u n \u n \ 


cf> 

x — n ? 


To understand the structure of this result, first notice 
that our problem conserves global charge. This is equiv¬ 
alent to the statement that it is invariant under a uniform 
static shift of all variables > <f>i + 0. Such shifts leave 
the fields $i 2 3 manifestly invariant, while the zero mode 
action for <J > 0 does not change provided 9 = const. 

Turning to the action of the remaining two fields $i i2 it 
is straightforward to verify that in the absence of Klein 
factors, i.e. for c-number valued scattering coefficients 
tij, the scattering Hamiltonian ([ 6 ]) becomes separable 
once Q = const, is frozen out. The action would then 


decouple into two independent actions of fields $ 72 , de¬ 
scribing forward ( 713 , 724 ) and backscattering ( 714 , 723 ), 
resp. These actions are obtained from S, Eq. ( |25| ), if one 
ignores the off-diagonal terms in the kernel of So- 

The Klein factors account for the non-commutativity 
of the scattering operators. In the reduced representa¬ 
tion, after integration over all auxiliary fields, their her¬ 
itage is the canonical contribution to the action (the 
off-diagonal terms in the matrix kernel) i.e. a term 
~ / dr $ i<9 T <i>2 stating non-commutativity of the forward 
and the backward scattering field. 

What are the commutation relations between the 
fields? Usually, the commutator between two conjugate 
variables [q,p] = ia is related to an action ^ f dtpq. 
Combining the two off-diagonal elements in the action 
Eq. (25), this would lead to [<J>i,$ 2 ] = mg. However, 
this result is incorrect. As usual in a functional integral 
based approach, operator commutation relations are ob¬ 
tained by analyzing the effects of time ordering in the 
integral, i.e. by subtracting the expressions ($i(±e)$ 2 ) 
from each other. The evaluation of these terms must take 
the presence of the diagonal dissipative contributions to 
the action into account. We then obtain, with e —> 0+, 

($i(e)$ 2 (0)) - ($i(-e)$ 2 (0)) 


r°° <hj_ 

— rv, 27 T 


'giwe _ g— icoe' 


2 U! 


ig / smw .irg 

— / auj -= * — 

2 i -00 w 2 


This implies the commutation relation 

]=i 7 f, 


(26) 


(27) 


which will play an important role throughout. 

Before ending this section, let us comment on a point 
that will play a role in the interpretation of our results 
below: the effective action (25) describes the physics of 


the so-called dissipative Hofstadter model introduced by 
Callan and Freed®] Recall that the Hofstadter model 
describes the physics of a rectangular lattice subject to a 
perpendicular magnetic field. Describing nearest neigh¬ 
bor hopping on the lattice in terms of two operators pi 2 
as cos(py 2 ), the action of the system contains a sum 
of two cos-operators, where different coupling constants 
may account for anisotropy. The presence of a mag¬ 
netic field implies the lack of commutativity, \p x ,Py] 7 ^ 0, 
i.e. a term ~ f drpid T p 2 in the action, whose cou¬ 
pling constant is a measure of the field strength. Fi¬ 
nally, dissipation may be introduced by including opera¬ 
tors ~ Y n \un\(\pi,n\ 2 + \P 2 ,n\ 2 )- Adding everything up, 
we arrive at an action equivalent to (25) (an identification 
Pi <-)• understood.) More specifically ( f25j ) describes 
the DHM fine tuned to a configuration where field and 
dissipation strength balance each other. What motivated 








Callan and Freed to generalize the non-dissipative Hof- 
stadter to the presence of dissipation was the expecta¬ 
tion that the fractal structure of its single particle spec¬ 
trum might turn into a fractal pattern of phase tran¬ 
sitions. Their analysis indeed confirmed that the phase 
plane spanned by field and dissipation strength is covered 
by a fractal network of phase transition lines separat¬ 
ing phases of localization (cos-coupling constants scaling 
to zero) from phases of de-localizatoin (diverging cos- 
amplitudes.) They arrived at that conclusion within a 
non-perturbative analysis resting on an approximate and 
an exact self-duality symmetry of the model. The above- 
mentioned balancing of field and dissipation strength in 
our context means that we are cutting through this pat¬ 
tern along a line parameterized by the Luttinger inter¬ 
action strength g. Later on, we will confirm within an 
instanton approach generalizing that of 3 ^ that for repul¬ 
sive interactions, g > 1, we are generically in a ‘delo¬ 
calized’ sector of the phase plane, which is another way 
of saying that the intermediate scattering fixed point is 
approached. 


from KF scaling with dimension 1 — <7 is that fluctua¬ 
tions in the incoming and outgoing chiral are partially 
correlated with each other. This is due to the fact that 
for energies below the charging energy of the dot, every 
incoming charge fluctuation is split equally between the 
two outgoing chirals, such that the amount of indepen¬ 
dent fluctuations is effectively halved. Thus, the pres¬ 
ence of dissipation makes the quadratic action unstable 
with respect to coherent scattering even for the hypothet¬ 
ical case of chiral edges with an ’’attractive” interaction 

1 < g < 2 . 

The analysis described above relates to the point ”1” 
in close vicinity of the Gaussian fixed point with = 0, 
7 „ = 0, denoted as (0,0) in Fig. 5. Due to the scal¬ 
ing dimension 1 — g/2 of coherent tunneling terms, the 
Gaussian fixed point (0,0) (WWFP) is unstable in the 
physical regime of g < 1. We want to emphasize that the 
variables y/j, are not conductances for tunneling in 
horizontal or vertical direction, but are instead the tun¬ 
neling amplitudes, which act as mass terms for the fields 
$1 and *f >2 in the regime of strong coupling. 


IV. THE WEAK-WEAK FIXED POINT: 

PERTURBATIVE ANALYSIS 

In this section, we will analyze how initially small co¬ 
herent tunneling terms Sh , S v will grow under the RG, 
depending on the value of the Luttinger parameter g. To 
leading order in the couplings 7 h and y„, only one of the 
fields $1 or $ 2 enters the analysis, and hence we only 
need the correlation functions 


• (28) 

2 \uj n \ 

This information allows us to work out the flow equations 
for jh, 7 „. Decomposing the fields <I>,; as usual into slow 
and fast modes and performing a RG transformation, one 
sees that the exponential of the field is renormalized by 
a factor 


e iV2<S>i b e iV2<S>i e -TY2 n ,f(®i,n$i,- n ) 


(29) 


Here, the factor b in front of the exponential is due to 
the rescaling of the integration measure dr in the action 
Eq. (25), and the J/f runs over fast frequencies A /6 < 
|w„| < A. Using the correlation function Eq. (28), we 
find 


I 1 

nj 




2 JA/b w 2 


(30) 


Taking everything together, we find that the scattering 
operators scale up with dimension 1 — g/2 , the result 
obtained earlier by us 11 . The reason for the deviation 


V. THE WEAK-STRONG FIXED POINT 
(WSFP): PHASE SLIP-INSTANTON 
MECHANISM 

We now want to assume that one of the two tun¬ 
nel couplings is much larger than the other one, say 
r = 7 „ » 7 = 7 ^. To heuristically understand what 
is happening in this limit, let us temporarily go back to 
a Hamiltonian description in which we have two tunnel¬ 
ing operators 

rcos(v / 2$ 2 ), 7C0s(v / 2$i), 

where the phase operators obey the commutation re¬ 
lations [$i,$ 2 ] = ign/2. To leading order in 7 , the 
interaction-picture time evolution of the exponential 
exp (—rcos(v / 2 < E > 2 )) under the action of the perturb¬ 
ing operator oc 7 is governed by the commutators of 
the type [ycos(\/2<i>i),exp (—r cos(v / 2<1> 2 ))]. The non¬ 
commutativity of $1 and $ 2 implies that the exponential 
of $1 acts as a translation operator for >I > 2 by an amount 
g-K. How does the system respond to such a transla¬ 
tion? A naive and incorrect argument suggests that such 
a phase flip by n (in the case g = 1 ) of \/ 2&2 will cost a 
large amount of energy due to the large magnitude of F, 
hence it should be accompanied by a second, almost im¬ 
mediate instanton that will add up — n or +7r to the first 
one, taking the field configuration of $2 back to a mini¬ 
mum. This argument suggests that the lowest terms ef¬ 
fectively contributing to the perturbation expansion will 
be of order y 2 . The less naive and correct heuristic argu¬ 
ment takes into account that the argument of cos(v / 2 <f> 2 ) 
contains finite k inodes besides the zero-modes generating 
the commutation relations of $ 2 relative to $ 1 . For this 
reason, the original phase change by gir due to the Klein 
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factors can be compensated by a kink in the finite mo¬ 
mentum part of $ 2 - This secondary kink enables the field 
$2 to settle back into an energetically favorable configu¬ 
ration without the help of another tunneling event, but it 
does cost a finite amount of dissipative action. The more 
quantitative analysis to be detailed momentarily shows 
that the primary effect of the added action contribution 
is a change of the scaling dimension of 7 to 1 — g, which 
should be compared to the dimension 1 — g/2 character¬ 
izing the weak-weak fixed point discussed in the previous 
section. 

On a more formal level, we will will apply an instan- 
ton analysis in <3> 2 ■ The complementary operator intro¬ 
duced within the duality transformation is assumed to be 
perturbatively weak. The instanton approach assumes 
that the relevant field configurations are rare phase slips 
(j> 2 = ^(2 1 + 1 ) 7 r —► ^(2(Z ± 1) + l)7r between 2 tt - 

consecutive minima of cos V 2&2 ■ Assuming that these 
events occur at times 77 , the phase profile is best de¬ 
scribed by its time derivative 

N 

d T $2 = ^2 V^ 7TS id( T ~ T i) > 

N 

^>n$2,n = ^2 'KS i e~ lU)nTi , (31) 

i =1 

where 5 is a broadened (5-function, and s, = ±1 denote 
the direction (the charge) of the phase slip. To facilitate 
the evaluation of the instanton functional integral, we 
use a Hubbard-Stratonovich with a field ©„ to partially 
decouple the quadratic term \co n \^ 2 ,n^ 2 ,-n according to 

T ^ ( (^ 27 j-g n || T ^^ < hi n $2,— ^ 

— t T f $1 n^l.-n\bJn | + —T-0n0-n|w„| 

„ \ n 9 2 

— n \^n | 0n^2,— n^n) • (32) 


Here, S 7 [$i] =7 f dr cos(v / 2 < ki), and 

S 0 [e,$i]=TV (—$!,„$!,_ n |w n | (34) 

V \ n 9 



We are now in a position to explore scaling dimensions. 
Due to the assumed weakness of both exp(— St) and 7 , 
we can treat the two operators independently of each 
other. The basic formula we use is that for a field 0 with 
correlation (</> n <^_ n ) = | , we have 

e ic<p 

in an RG step. In our particular case, (X n X_ n ) = 
yjT—j-M -1 , where M is the matrix controlling the free 
action: 

($l,n$l ,-n) = rr,} T TTg => K = 7 Tg,C=V2 

1 | id n | 

(0„0_„) = 1 . — => k = —, c = a/ 27r. (35) 

T\u n \ 7 xg 7 xg 

We then find that the two coupling constants scale as 

e~ Si —> e _Si 6 1 _ 9 , 7 —> 7 b 1-9 . (36) 

Thus, in the regime g < 1, the weak coupling 7 is a rel¬ 
evant perturbation, while the exponential of the instant 
on action scales to zero, indicating that instantons in the 
strongly coupled scattering channel are irrelevant. The 
weak-strong FP has an unstable direction corresponding 
to the growth of the "weak” scattering channel, while 
the irrelevance of instantons in the ’’strong” channel in¬ 
dicates that the channel stays at strong coupling. This 
situation corresponds to point ” 2” (and similarly to point 
”4”) in the phase diagram Fig. 6 . 

VI. THE STRONG-STRONG FIXED POINT 
(SSFP): INSTANTON ANALYSIS 


Doing the Gaussian integral over 0, one gets from the 
second to the first line. We note that 0 and $2 are a 
canonical pair. Denoting the instanton action by S), the 
functional integral now becomes 


Z = j D$iD 0 e " So[e ’ <I,l] " 5 ^ [ ' i>l1 


(33) 


X E^lE /IW 7VS<e ^ 7 r Ef=l5 ’ TE " 

N '{«<} 


e„e“ 


J D^iDQe~ S °^ & ’^ +e Si / dTCOS (v / 2 u-e)-S 7 [$i] 


We now want to discuss the situation when both the 
vertical and horizontal tunneling terms are strong, such 
that none can be discussed perturbatively anymore. In 
the ground state of the system both fields and $2 
then occupy one of the minima of the cos-terms, and 
elementary excitations above the ground state are jumps 
from one such minimum to the neighboring minimum, 
see Fig. 7. In this way, the field configurations move 
on a lattice of ground states, which are connected by 
instantons (jumps) in the fields. The scaling dimension 
for such a jump can be obtained by computing the action 
for the instanton due to the quadratic dissipative term in 
the Helds, as described in the section above. In this way, 
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0 2tt $1 


FIG. 7. Square lattice for instantons in the (4>i, 4>2)-fields. 
Jumps in the horizontal and vertical direction have equal 
length, and instantons in diagonal direction are less relevant 
than horizontal and vertical ones. 


potentials. Assuming these events to occur at times 7 y a , 
the phase profiles are best described by their time deriva¬ 
tive: 


Ni 


d T $i = ^2 27rS i,afi( T - n,a), 

a—1 

Ni 

iu n §i, n = ^ 27 TSi^ a e~ lUnTi ' a . 


(37) 


where S is a broadened <5-function, and Si ja = ±1 denote 
the direction (the charge) of the phase slip. To facili¬ 
tate the evaluation of the instanton functional integral 
we Hubbard-Stratonovich decouple the quadratic depen¬ 
dence in d T 4 »j or iu> n $i in . Writing the free Gaussian ac¬ 
tion as S'o[ < f ) ] = j J2 n 1 ^-n> where the matrix 


X is defined through Eq. (25) as X n = tt g 


( sgn(n 
= 7TO . 


V -1 sgn(n 


We complete the square in terms of a dual integration 
variable 0 as 


we establish that the strong-strong fixed point (SSFP, 
denoted by the point 3 in Fig. 5) is indeed stable. In the 
following, we will present a more formal derivation of this 
result. 

If both > 0 are strong, we may subject the prob¬ 
lem to an instanton duality mapping. In this regime, the 
relevant excitations are rare phase slips ( F, = ( 2 Zj + l) 7 r —► 
( 2 (li ± 1 ) + l) 7 r between consecutive minima of the cos- 


So [©,$] = 

= ^2 ^n^>nX n Q- n + T Q^OJ n ^- n . (38) 

n n 

The second term tells us that $ and 0 are canonical 
variables. 

We denote the instanton actions as S\_ t and, keeping 
in mind their exponential smallness in the coupling con¬ 
stants, 'y Vt hi write down a double expansion 


2= //»«- tE - e -"- X - 9 -EwW7 F f 

Ni 1- 2 '{s i , a } v ' 


= / DQe 2 ^' 


- ^2 n QrL^nXnQ - n —'22 i C f d,T COS (©i) 


(39) 


r 


In the last step, we scaled 0 by a factor 7 r for convenience. 
Substituting the form of the X-matrix, we arrive at an 
action 


s[e]-££E«S(f: 


I 

n \^n 


Q-r 


(40) 


+ ^e Si ’* I dr cos( 0 i). 


dual to 
stants: g 


(25). The differences are in the coupling con- 
2 / 5 , and r ) v /h e~ Si,2//1 . This means that 
for very large and nearly equal coupling constants 7 j, the 
coupling constants of the dual theory can be renormal¬ 
ized by first order RG analysis where the reasoning of the 


foregoing section obtains the scaling 

e -SM_> e -SM 6 l-*. (41) 

At first sight, this result looks surprising: for g < 1 the 
instanton insertions are irrelevant, and the system flows 
to strong coupling, while for 1 < g < 2 the instanton op¬ 
erators are relevant, and the weak potential perturbation 
addressed in the previous section is also relevant. This 
suggest the existence of an intermediate fixed point. Fi¬ 
nally, for 5 = 1 , the instanton operators are marginal, a 
situation which seems inconclusive at first. 

In order to clarify this situation in more depth, we 
follow Callan & Freed® and generalize the matrix X n as 


x -i _^ J_ ( a sgn(n) 


_ 13 

2i t V -/? a sgn(n) 


( 42 ) 
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i.e. we allow for independent coefficients of dissipation 
(a) and canonical term (/3). Our current situation is 
recovered by setting a = /3 = K With this generaliza¬ 
tion, we have ($»,„$*,-n) = and the P er ~ 

turbative scaling dimension becomes 1 — • This 

tells us that the weak scattering operators are relevant if 
a < a 2 +/3 2 . Defining z = a + i/3, we realize that z must 
lie outside a circle in the complex plane centered around 
(1/2,0) and with radius 1/2 (cf. Fig. [8}) In our case, we 
are sitting on a line z = P(1 + i), and for g < 2 we are 
outside that circle, i.e. we are dealing with relevant op¬ 
erators. We note that for /3 = 1, the condition collapses 
to g < 1, corresponding to the standard KF case. 

In the dual case, we have the same situation, only that 
X- 1 gets replaced by X" 1 -7 ^ sg n(n) 

Defining z = a2 (/«2 — i a 2 +a 2 > we realize that the instan- 
ton duality amounts to a mapping z 5 = z -1 . The 
instanton operators are irrelevant, provided z lies inside 
the circle above. In our present setting, this requires that 
l/(s _1 (l + i) = f (1 — i) lie inside the circle. For g < 1 
this is the case. In the hypothetical case of attractive lead 
interactions, g > 1 we are outside and the instanton op¬ 
erators are relevant, while the weak scattering operators 
are also relevant. 

For g = 1, we are borderline. We note that backscat- 
tering is a relevant perturbation in the weak-weak limit 
- with a scaling dimension of the backscattering proba¬ 
bility equal one. This is in line with reference [311 (cf. 
Eq. (346) with N c h = 2, see also reference hod- In the 
strong-weak and in the strong-strong limit the flow of 
the backscattering probability is marginal. This is in¬ 
deed verified in the next section (cf. also Fig. 10). For 
this behavior to hold, we have utilized in our analysis the 
instantaneous character of backscattering. Reference [31] 
allows for more general scenarios which may lead to trun¬ 
cation of the scaling flows (cf. Eq. (348) there). 

In the next section, we treat the model for the spe¬ 
cific value g = 1 of the Luttinger parameter by us¬ 
ing the method of fermionization, and find agreement 
with the scenario discussed here. Even if we consider 
non-interacting leads, yet, there is still interaction in 
the model: the electrostatic charging energy, which cou¬ 
ples the four chirals among each other. It is clear from 
the scaling obtained above that the flows in the vicinity 
of both the WSFP/SWFP and the SSFP are marginal. 
This nicely reproduces the expectation that for Fermi- 
liquid leads, scattering should be marginal in the low- 
energy limit, and that arbitrary ratios between forward 
and backward scattering strength can be realized. 

VII. GENERAL COMMENTS ON THE 
ANALYSIS 

We would like to discuss some subtleties of the analy¬ 
sis presented above: (i) We assume a cutoff (voltage or 
temperature), which may be small, yet larger than the 



FIG. 8. Relevancy of perturbations as a function of the two 
parameters a,/3 defined in the text. The white disk/shaded 
area is the region of the phase diagram where scattering terms 
are irrelevant/relevant. At the weak coupling fixed point, a = 
/3 = g- 1 and the region of relevant perturbation (shaded) is 
reached as g \ 2 (line departing from origin corresponding to 
lowering g). At strong coupling a = —/3 = g /2 and the region 
of irrelevant perturbation (circle) is reached for g \ 1 (on the 
line /3 = —a.) The line /3 = 0 corresponds to the KF model, 
where perturbations around weak coupling become relevant 
at g = 1. This means that for repulsive lead interactions, 
g < 1, perturbations at weak (strong) coupling are relevant 
(irrelevant) implying the stability of the stability of the strong 
coupling fixed point. 


level spacing A of the QD (which for a typical quan¬ 
tum Hall geometry may be exceedingly small.) Had we 
pursued the RG all the way down to A, the character 
of the model would have changed: there are no inelas¬ 
tic excitations then, hence the inelastic channel is frozen. 
In that case, in the presence of only elastic scattering, 
we are driven to either full transmission or full reflec¬ 
tion, in agreement with the Kane-Fisher result. In the 
present analysis, we formally take the limit of an infi¬ 
nite QD, before pursuing our RG procedure all the way 
down to zero temperature or zero voltage, (ii) In the ab¬ 
sence of forward scattering channels (714 = 723 = 0 ), our 
model and the results obtained for its scaling reduce to 
those considered by Furusaki-Matveev. (iii) Except for 
the highly degenerate points of perfect scattering ampli¬ 
tude cancellation 7 13 + e _l 7 r 9 7 2 4 = 0 or 723 + 714 = 0 , 
our analysis holds regardless of scattering phase shifts 
and strengths (at the degenerate point, the model col¬ 
lapses to a variant of the Kane-Fisher problem). How¬ 
ever, these points require fine tuning and are nongeneric. 
The specific conditions for perfect scattering amplitude 
cancellations are related to the choice of the Aharonov- 
Bohm phase enclosed by the the corresponding tunnel¬ 
ing path. Here we have followed special choice, that of 
a zero Aharonov-Bohm flux. We note that OCA Ref. 
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(in their 3 lead geometry) have allowed the freedom of 
choosing the Aharonov-Bohm flux, (iv) An extension of 
our analysis from FQHE to interacting Luttinger wires is 
presented in Appendix 


VIII. REFERMIONIZATION FOR THE CASE 

9= 1 

We consider a special model, where there are only two 
nonzero scattering amplitudes 714 = 7 h and 724 = 7 „, 
see Fig. [9j and where the leads are non-interacting. The 
model has an intereting dynamics all the same, since at 
low energies the mode $3 describing charge fluctuations 
in the quantum dot is already frozen. In this limit, the 
tunneling amplitudes 7 h and 7 „ are both relevant pertur¬ 
bations. Using a mapping onto non-interacting fermions, 
we derive an exact solution of the model, which describes 
the competition between the competing scattering pro¬ 
cesses 7 h and 7 „. We find that at high temperatures and 
small 7 h and 7 „, both 7 /j and 7 „ increase tx T -1 / 2 , in 
agreement with a weak coupling perturbative RG in Sec¬ 
tion IV. The asymptotic value for T —> 0 is non-universal 
and depends on the ratio of initial values of 7 h and 7 „. 
This finding complements the results of Sections [V] and 
|VI[ where it was found that non-interacting leads with 
g = 1 separate two different regimes of RG flow. 



1 

\ 
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FIG. 9. Setup with two competing scattering amplitudes 
714 = 7 h and 724 = 7 u- We consider the special case of non¬ 
interacting leads g — 1 , and denote the probability for elastic 
scattering from segment 1 to segment 4 by P 14 , and the prob¬ 
ability for elastic scattering from segment 2 to segment 4 by 
P 24 . 


The starting point for our discussion are the tunnelling 
terms and the dissipative term in Eq. (25). Due to the 
off-diagonal terms in So , the fields $1 and *I >2 satisfy the 
commutation relation ['f > i, < f > 2 ] = 77r|. In this appendix, 
we specialise to non-interacting leads with g = 1. Due to 
the presence of the charging energy on the dot, the diago¬ 
nal entries of Sq describing dissipative processes are given 
by |w n |/27r. We compare this value with the dissipa¬ 
tive action |w n |/ 47 r^ for the charge density field describ¬ 


ing backscattering of quasi-particles (in the absence of a 
QD) between two counter-propagating fractional quan¬ 
tum Hall edges with Luttinger parameter v. We find 
that non-interacting leads with a QD charging energy 
have a dissipative action equivalent to that of interacting 
v = 1/2 LL edges in the absence of a QD. The special case 
of backscattering between v = 1/2 LL edge states is can 
be solved exactly by the method of refermionization 3 ^, 
which we will adapt to our problem of two competing 
scattering processes in the following. In order to achieve 
this goal, we undo the shift Eq. (B.17), such that in the 
following the fields <f>i and < I > 2 contain finite frequency 
modes only. The non-commutativity between the fields 
is now described by Majorana operators 771 , 773 , and 774 
with {gi, ? 7 j} = 2Sij. Switching to a Hamiltonian formal¬ 
ism suitable for refermionization, the tunnelling Hamil¬ 
tonian is given by 


e iV2<S>i(0) + e -iV 24>i(0) 


H tun = immih 

+77727747, + e -^2(°) 


(43) 


Here, we have absorbed a factor 1/2 into the tunnelling 
amplitudes 7 h and 7 v as compared to Eq. (25). Tak¬ 


ing into account that the dissipative actions of $ 1 , $2 
are equal to that of an infinite chiral Luttinger liquid, 
the operators e *^ 2 ® 1 / 2 can be interpreted as the opera¬ 
tor fI/ 1/2 of free chiral electrons. However, to make this 
equivalence more precise, we need to introduce Majorana 
fermions fi/ 2 as a Klein factor for the new electron^!!, 
and define 


e n/2$ 1/2 (0) = ^ 2 ™* 1/2 (0)/ 1/2 . (44) 

Here, a denotes the short distance cutoff of the theory, 
which is needed to make sure that the fermion fields ^ 1/2 
have the proper dimension of one over square root of 
length. In this way, we obtain the fermionized tunnel 
Hamiltonian 


H t un = 777 2 7747 1) V27ra 


* 2 (0) h + h *5(0) (45) 


+ir]imihV2TTa 'Li(O) /1 + /1 ^1(0) 


In addition, the fermions have a kinetic term (the velocity 
is taken to be unity) 


Hr, = dx 


*l(z) (-id x ) *i(a;) + 4)\(x) (-id x )4> 2 (x) 


(46) 

Since the leads are non-interacting, we can discuss 
transport through the QD in terms of a scattering for¬ 
malism. For an incoming particle on edge 1, we define the 
probability P 14 to scatter onto the outgoing edge mode 4. 
Similarly, for an incoming particle on edge 2, we denote 
the probability to scatter onto the outgoing edge mode 4 
by P 24 . Then, in Appendix D we derive the result 
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FIG. 10. Description of the RG flow in terms of scattering 
probabilities Pi4 and P24, which in the weak coupling limit 
are proportional to the square of tunneling amplitudes 7 ^ and 
7 ^, respectively. Near the RG-unstable weak-weak fixed point 
(red circle), there are perturbative corrections to the Gaussian 
(1/2,1/2) conductance in both 7 ^ and 7 „, which grow under 
the RG. The RG flow terminates on a line of fixed points 
(red line), which are characterized by scattering probabilities 
and a conductance which depend on the ratio 'iv/'lh of bare 
scattering amplitudes. 

Pu = y f ( f ) ( 47a ) 

p ™ = Y F (j) ’ (47b) 

with 

T h = 47T07^ , T v = Airajl , T = T h + r„ , (48) 
and 

2 oc 

F(x) = - arctan — . (49) 

The scaling function has the limiting behaviors F(x) = 
1 — |(*/2) 2 for small x, and F(x) = 2 for large x. Using 
these asymptotics, we find 

( % h for T > T 

Pu = \ , ( 50 ) 

[ ^ for £ -> 0 

and similarly for P 24 , with Fft and F„ interchanged. 
In the high temperature limit, this result agrees with 
the perturbative analysis presented in Section |TV| Since 
Pi 4 oc 7 ^ in the high temperature limit, the flow of P 14 
is representativ for the flow of the tunneling amplitude 
7 /,, and similarly for P 24 and j v . The flow of Pi 4 and 
P 24 is shown in Fig. [lO] For weak bare values of 7 ft, 7 „ 
with Tft/T <C 1, I\,/T -C 1, the flow starts in the vicinity 


of the unstable weak-weak fixed point. In terms of the 
variables P14 and P24, the flow occurs along straight line 
trajectories, which stop at the line of fixed points defined 
by P14 + P24 = 1. The end point of the flow on this fixed 
line is determined by the ratio of initial couplings jh/jv 
In particular, in the weak-strong limit with, say, j v 7 ft, 
the flow is truncated at a temperature T ~ r. u , and the 
weak coupling 7 ft is marginal. This situation is interme¬ 
diate between the case g > 1 , where the weak coupling 
is irrelevant, and the case g < 1 , where the weak cou¬ 
pling is relevant. Thus, the exact solution for the case 
of noninteracting leads with g = 1 is in full agreement 
with the perturbative analysis presented in Section [V] 
We find that for g = 1 there is no strong-strong fixed 
point, which again is in agreement with the analysis in 
Section VI in which it was found that g = 1 separates 
the case of a stable strong-strong fixed point for repulsive 
lead interactions with g < 1 from the case of an unstable 
strong-strong fixed point for g > 1 . 


IX. SUMMARY 


In summary, we have considered a model of scattering 
between one-dimensional chiral leads in the presence of 
both inelastic and elastic scattering channels. We have 
found that this combination of scattering can stabilize 
a new \ - \ fixed point, with probability 1/2 for the 
transmission and probability 1/2 for the reflection of an 
incoming particle. This \ \ fixed point is intermediate 

between the 1-0 and 0-1 fixed points in the presence 
of only elastic scattering channels. In order to estab¬ 
lish the existence of this fixed point, we employed a non- 
perturbative instanton analysis in either one or both of 
the scattering channels. Our main result is that the in¬ 
termediate fixed point is stable for Luttinger parameters 
g < 1 in the leads. This conclusion is backed up by 
a refermionization analysis for the special value g = 1 . 
For non-interacting fermions with g = 1, we recover the 
well-known marginal relevance of scattering. While this 
1/2 — 1/2 fixed point is the main result of the present 
analysis, we recall that our previous paper 11 pointed out 
novel results concerning noise and current-currnet corre¬ 
lations. 
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Appendix A: Proof of Eq. (161 


As usual, the path integral fixes operator ordering 
through time order. Consider, for example, two quasi¬ 
particle operators, 

4q(r + <5r)4' i (r) -a e^ 2 *' 
xe *(f 

x^i(r + . (A.l) 

The expression we want to compare with looks the same, 
except for an exchange of time arguments (r — St) <->■ r. 
Now, imagine this expression inserted in the functional 
integral. We may always integrate over {qi}, as these 
fields enter the action linearly (assuming a perturbative 
approach, where all dq. are expanded out of the expo¬ 
nent.) Integration over the g’s generates the step func¬ 
tion profiles 


Pk (T) =Pk,0+ g^SqQjT -Tg), 
a 

where pkp is the initial value, the sum runs over all ap¬ 
pearances of qk in the exponents, r a are the respective 
times, and s a = 1 for a — iq and —1 for a +iq factor. 
Now, with these structures in place, we can explore the 
behavior of the exponent above. Denoting by <3> ? y the 
cumulative phase which will not respond to an exchange 
of the time arguments, we have 

^i(r + 5 r)$j(r) -A , 

*i(r)Vj(T + St) -a 

where the phase in the first line reflects the fact that at 
time r pj jumped by g/V 2 and this phase change is read 
out by the contribution a.ijPj{r + St) to the first phase. 
Similarly for the second line. The relative phase is given 
by 

— a ji) = ^)' 7T 9 a ij 


as required. 


Appendix B: Derivation of the effective action (251 


In this section, we derive the effective three-variable ac¬ 


tion (25) from the original description in terms of twelve 


fields, (j>i,qi,Pi . We start by casting the orthogonal trans¬ 
formation (23) into the more compact notation 


/I 1 1 1 

*= M =n I:! -i "i 1 

Vi l -l -1/ 


This transformation is orthogonal, M T = M 1 , and thus 
the inverse transformation is given by <f> = M _1 $ = 


M t $. For the zero modes, this allows to define p = 
M~ 1 P and q = M~ X Q as above to obtain 

/n 

Q = -2 P 0 - —$3- (B.l) 

7r 

Then, the zero mode part of the scattering operators as¬ 
sumes the form 

(-D+SPol^ig-ifs, 

§ 4 4-2 = e K+ Q 2-Q 3 ) e i^ f+D-Pol^^gifs, 

4 > 3 \I >2 = e *(+Qi~Q 3 ) e i% (- p 2+ p o)^ 3 ^ 2 e“^ 9 , 

'F 4 'F 1 = gH-Qi-Qalgif (+ p 2+Po)^ 4l p ie -i§9_ (B.2) 

The transformations (q,p) -A (MQ, MP) is orthogonal 
which means that the canonical piece of the action re¬ 
mains unaltered, as 


^ J dTPid T Qi. (B.3) 

In principle, there is also the zero mode charging term 
H ~ J2i Pf to consider, however, it does not play any 
role due to the prefactor of inverse contour length, which 
becomes zero in the thermodynamic limit. 

We note that the scattering terms do not couple to 
Qq. This means that we may do the integral over Qq 
to conclude that d t Po = const.: the total charge on the 
system is constant. Without loss of generality, we may 
call this constant 0. Our scattering terms thus simplify 
to 

'F 3 'F i = e*(-Q 2 -Q 3) e -if- Pi ^3^i e -ifs ) 

ij/ 4 «F 2 = gil+Qs-Qalg+if 

4-3^2 = ed+Qi-Qalg-if ^s, 

§ 4 Ti = eil-Oi-Oslg+if Pj^ig-if (B.4) 

Next, ignoring the overall charging term, we observe that 
P 3 does not appear in the scattering operators, nor in the 
dot charging term. This leads to the constraint d t Qj, = 0, 
and we may set Q 3 = 0 without loss of generality. As a 
result we have the further simplification, down to 


T 3 Ti = ig-ifs, (B.5a) 

vj/ 4 <F 2 = e +iQ2 e +i ^ Pl ^ 2 e^ 9 , (B.5b) 

T 3 T 2 = (B.5c) 

# 4^1 = e~ iQl e +i ^ P2 ^ie-^ 9 . (B.5d) 


We note that the zer o mode sectors of t he sca tt ering oper- 
ators in Eqs. (B.5a|, (B.5bl and Eqs. (B.5c), (B.5dl now 
pairwise commute among themselves. However, those in 
Eqs. (B.5a), (B.5bl and Eqs. (B.5c), (B.5dl do not com¬ 
mute. This structure suggests a further simplification: 
from the commutation relation 


[Qi, Pj] = giSij (B. 6 ) 
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we compute 

[Q 2 + |-Pi, Qi - |-P 2 ] = -ign. (B.7) 
Now this relation motivates the canonical transformation 


A = -^ + ^Pi, 

-v/'TT 2 


x , = ft + 0 


2 , 



+ 



2 


P 2 


y' = 




(B. 8 ) 


The new variables obey the relations [ X , Y] = ig and the 
same for X ', Y'. The primed and unprimed variables are 
mutually commutative. All this means that the canonical 
piece of the action remains invariant. In particular, we 
have a contribution 


S[A,F] = l - J drYdtX . 

The scattering operators assume the form 

T 4 T 2 = e i 'f* x, ij)4'ip2e i ^ 9 , 
T 3 T 2 = 

T 4 T x = e i v' 5Y $ 4 V> 1 e- < 5® ) 


(B.9) 


(B.10) 


which means that we do not need the variables X',Y', 
and do not consider them in the following. 

Using the canonical representation (23) and assuming 
locking of the charge mode we have 


^ 3^1 nj e i(-YnX+Y2$ 2 ), 

$3^r 2 ~ e H-V* Y -V2*i) t (B.ll) 


where we have absorbed the exp(±iv/ 7 r/ 2 <?) factors in the 
scattering phases 713 and 723 . 

The structure of the scattering operators suggests a 
shift 


$2 -A $2 + \j\ X i 

*!-►*!- \J\y, (B.12) 

whereupon the scattering operators simplify to 

vp 3 vE - 2 ~ e _iv/ 5 $ i, (B.13) 

As a result of this shift, the fields *I>i and <f > 2 no longer 
commute with each other, and we have 

[$ 2 ,$i ] = w\ ■ (B.14) 

Now we are in a position to express the local effective 
action in terms of variables which allow an efficient analy¬ 
sis of the competing scattering processes. The dissipative 


part of the action, obtained by integrating out the finite 
wave vector modes of the bosonic fields, is given by 


Sdiss[$] = ^ EEi 

^ i —0 n 


1$,; 


(B.15) 


where d>i ira denotes the n-th Matsubara component of the 
field The sum runs only up to i = 2, on account of 
the locking of the held <J> 3 . The i = 0 contribution to 
the sum defines the zero mode action (24), and will be 
ignored throughout. Introducing the shorthand notation 


$ = (<f-i,$ 2 ) 5 


= (X,F) 5 


(B.16) 


the quadratic part of the action then assumes the form 


So [X 4>] = — £ - -- E $ 

ng z ' g n z ' 


t — I 

n\^n 


T 


/ |u 

2 g ^ \ ~ 


UJ n 


\ ^ 


W n \Un 


(B.17) 


The coupling between 4> and H makes sure that the com¬ 
mutation relation Eq. (B.14) is reproduced in correlation 
functions of the <£ 7 . Since the vector variable H no longer 
appears in the scattering part of the action, it can be 
integrated out. Doing the Gaussian integral, we obtain 
the action So[$i,$ 2 ] as given in Eq. (25). 


We finally note that the two ’vertical’ scattering oper¬ 
ators ~ 4 f 3 ^'i/'I' 44'2 (and similarly for the horizontals) 
may couple at arbitrary strength/scattering phase. This 
leads us to 5[3>] = S'o[$ 1 ,$ 2 ] + S v [$ 2 ] + <S)J$i], where 
the scattering action reads (cf. Eq. (fl9[)) 


Sv[$2]= Y, [ dr | 7 X | cos(v / 2 < I >2 + </> x ), 

x= 13 , 24 "' 

Sh[®i] = E [ l 7 x| cos(v / 2 < I > 1 + </> x ). (B.18) 

x= 14,23 •* 

Here, the phases <pi : absorb the phase of the complex 
scatteri ng am plitudes 7y, the phases exp(* 7 T< 7 / 2 ) appear¬ 
ing in (B.17), and the relative sign change of the field 
variables in 13 vs 24 (cf. Eq. ( B.10[ )). 

Without loss of generality, we can use an addition 
theorem for trigonometric functions followed by a con¬ 
stant shift of the fields ‘Iq and <h 2 1 ° transform the 


scattering actions into the form given in Eq. (25), with 
7v = | 7 i 3 + e _I7r9 724I and j h = (723 + 7 i 4 |- 


Appendix C: Non-Chiral Luttinger Liquids 

Throughout the analysis presented in this paper, we 
have assumed a setup consisting of a quantum dot con¬ 
nected to 4 chiral wires, which are geometrically sepa¬ 
rated from each other. These chiral wires talk to each 
other through the tunneling to the QD (and the charg¬ 
ing interaction thereon), and through elastic tunneling 
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terms between the incoming and the outgoing channels. 
Such channels can be realized, for example, as the edges 
of a fractional quantum Hall strip; the most relevant tun¬ 
neling operators are those of fractionally charged quasi¬ 
particles (the latter possess fractional statistics as well), 
and we have assumed that such tunneling terms are al¬ 
lowed. While intra-(chiral)channel interaction is allowed, 
inter-channels interaction is excluded. 

When it comes to realizing our theory with non-chiral 
Luttinger liquid wires, the situation is trickier. There are 
two main issues that should be noted. First, through¬ 
out our analysis, we have assumed that forward and 
backward (elastic) scattering are treated on equal foot¬ 
ing. There is no concrete significance associated with 
forward or backward. When it comes to non-chiral Lut¬ 
tinger wires, unless special conditions are specified, for¬ 
ward (elastic) scattering is irrelevant, and is clearly dis¬ 
tinct from backward scattering. Second, as was noted 
above, in the case of chiral edges supported by an incom¬ 
pressible fractional quantum Hall electron gas, it is clear 
that tunneling of quasi-particles is allowed. The situation 
with non-chiral Luttinger wires is different, however. In 
the analysis depicted below we allow only for electron 
tunneling; in general the results of the ensuing analy¬ 
sis will be qualitatively different. However, for a specific 
choice of the interaction (compare Eq. ( C.10| )), our beam 
splitter can be realized using non-chiral Luttinger liquid 
wires. 

Let us briefly review these two types of processes. We 
briefly repeat the analysis of Ref. [28] . Consider the Lut¬ 
tinger liquid Hamiltonian 


h = —[,K{d x (j>f + K~\d x ef] 


(C.l) 


[vr{x),Vr{ x ')\ = -[tp L (x),ip L (x')\ = iirK sgn(x-x'). 
The operators e lv>R / L are field operators of chiral quasi¬ 
particles that carry charge Ke. However, experimentally 
realizable tunneling processes can only involve elctrons, 
and in the following we will discuss the RG relevance of 
such processes. 


1. Interactions within one LL only 


A sketch of a possible realization is shown in Fig. 11 


It is well known that that a local potential (or, in the 
microscopic model, a weak bond) gives rise to backscat- 
tering and is a relevant perturbation in the RG sense. 
More problematic is the microscopic realization of an 
RG-relevant forward scattering term, which coherently 
transfers electrons across the quantum dot (dashed red 
lines in Fig. 1). In order to compute the scaling dimen¬ 
sion of such a forward scattering term, we describe the 
LL by chiral bosonic eigenstates ip± with imaginary time 
action 


Sl/R= 4jtK J dxd,T ( ±i9 r + d x) <Pr/l ■ ( C - 5 ) 

Here, u denotes the renormalized velocity, K the LL pa¬ 
rameter, and the smooth part of the total electron density 
is given by p(x) = ^(d x tp R + d x tp L ). In this basis, right 
and left moving electrons (i.e. electronic states at the left 
and right Fermi point of the noninteracting system) have 
the bosonized form 

V>r/l(z) = — ^ e i[K±VH{x)-K TVL {x)] s ( C . 6 ) 

y/2ira 


where the bosonic fields satisfy a Kac-Moody algebra 

[cj){x),e{x')\ = ySgn^-a/) . (C.2) 

Here, K is the Luttinger liquid interaction parame¬ 
ter. One can define left and right moving chiral fields, 
®r/l = 4> ± 9, and respective electronic held opera¬ 
tors ipR/L = e izlkFX e l ® R / L , with commutation relations 
[$ r ( o ;), $«(£')] = -[$ L (x),$ L (x')] = insgn(x - x'). 
In terms of these chrial modes, the Hamiltonian reads 

H = ttvq [p 2 R + p 2 L + 2A p R pL\ , (C.3) 

with an inter-channel interaction (between the two chiral 
modes). Here, 

Pr/l(x) = ±2-d x ^ R / L . (C.4) 

One can define other modes, <Pr/l = K(j)± 9 1 in terms of 
which the Hamiltonian decouples into ’’left” and ’’right” 
sectors. These new modes have commutation relations 


with K± = (K 1 ± l)/2. Thus, a local backscattering 
operator is given by 

0(x i) = (C.7) 

_ ^ c i(K+-K-)\<fi R (xi)+<PL(xi)] 

2na 

Using K + — I\ l = 1, we reproduce the well known result 
that the correlation function (Oback(i)Oback(O)) ~ t ~ 2K , 
which makes backscattering a relevant perturbation for 
a repulsive interaction with K < 1. On the other hand, 
a forward scattering term involves electron operators at 
two different spatial positions 

df(xi,x 2 ) = V’rO^V’rOo) (c.8) 

2 ira 

Since positions X\ and x 2 are separated by an ’’infinitely 
long” dot region, correlation functions between fields at 
positions X\ and x 2 vanish. As a consequence, the corre¬ 
lation function in time of the forward scattering operator 
factorizes according to 
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(Of(t)Of (0)) (X e -h K +^R.{xi,t)-Vn.(xifl)¥)-hK 2 Jr {[‘P'R.{x2,t)-‘PR{x2fl)]' 1 ) 

a f -KT(2K=+2^) 

« H K " 1+i ^ . (C.9) 



FIG. 11. It is not possible to realize a beam splitter with 
two standard non-chiral LLs because backscattering between 
left and right moving electrons in the same LL (short blue 
lines) is relevant, whereas forward scattering across the dot 
(dashed red lines) is irrelevant. The part of the wires inside 
the dashed box constitutes the quantum dot with a nonlocal 
charging interaction. 


We obtain V\ = diag(l + 417,1,1,1). In the following, we 
introduce the charge mode velocity 

u = 1 + 4U . (C.12) 


Next, we rescale fields such that the velocity matrix be¬ 
comes equal to the unit matrix, i.e. V 2 = AViA with 
A = diag(l/v/u> 1,1,1). With the help of these transfor¬ 
mations, we find K 2 = AAKA T A with 


/ 0 l/y/2u l/y/2u 0 \ 

l/y/2U 0 0 l/y/2 

l/y/2u 0 0 -l/y/2 

\ 0 l/y/2 — l/y/2 0 


(CAS) 


Since the function f(K ) = K _1 + K has its mini¬ 
mum value 2 at K = 1, the forward scattering operator 
Eq. (C.8) seems is irrelevant for all values of K. As a 
consequence, finding a microscopic description of an RG 
relevant forward scattering operator requires a modifica¬ 
tion of the interaction term. 


Since K 2 is real symmetric, it can be diagonalized by an 
orthogal transformation according to Kd = BK 2 B T with 


/ - I / a /2 1/2 1/2 0 \ 

l / y /2 1/2 1/2 0 

0 - 1/2 1/2 l / y /2 

V 0 1/2 -1/2 l/y/2 j 


(C.14) 


2. Local interactions between both nonchiral 
Luttinger liquids 


Here, we want to allow for the possibility that there 
is local interaction term including both LL wires as 
shown in Fig. [l2j Following the notation of Wen 
(Quantum Field Theory of Many-Body systems, sec¬ 
tion 7.4.6), we describe the system by a JF-matrix with 
K = diag(l, —1,1, —1), where we imagine that the first 
two branches belong to the first LL with fields (fiiR^o and 
ipiL, o, and the second two branches belonging to the sec¬ 
ond LL. These are eigenmodes in the absence of inter¬ 
actions. A local charging term acting on both LL wires, 
i.e. equally on all four modes is described by the velocity 
matrix 


V = 


1 + U U U U \ 

U 1 + U U U I 
U U 1+U u 
U U U 1 + U J 


(C.10) 


and we obtain Kd = diag(—1 /y/u, 1 /y/u, —1,1). The en¬ 
tries of Kd are the inverse of the scaling dimensions of 
the new fields. Through this series of transformations, 
the new fields 4>i , (j) 2 , </>3 , </>4 are related to the original 
fields according to the transformation (j) = BA~ l Atp. In 
order to compute scaling dimensions of tunneling opera¬ 
tors, we need the inverse transformation which is given 
by 


C = (BAA)- 1 = A t AB t = 


W \ 


— 1+ \/u 

l+Vu 

0 

1 

2y/2u 

2 \+2u 

n/2 

1 + +U 

2y/2 u 

1 — y/u 
2y/2u 

1 

V2 

0 

1+x/ii 

2y/2u 

1 -\-\^LL 
2+2u 

0 

1 

%/2 

l-\-y/u 

1 — y/u 

1 

0 

2\/ 2 u 

2y/ 2u 

V2 


(C.15) 

Using the matrix C, we can express the original fields in 
terms of the new ones according to 


We now diagonalize the real symmetric velocity matrix 
by an orthogonal transformation according to V —> V = 
AVA T with 


/ 1/2 1/2 1/2 1/2 \ 

I/a /2 -1/V2 0 0 

0 0 1/V2 -1/V2 

\ 1/2 1/2 - 1/2 - 1/2 J 


<p = C(f> 


(C.16) 


In order to compute the scaling dimension of the operator 
for backscattering a right mover in wire one into a left 
mover in wire two, we use the expression 


-rh t 


1 


03 + </>4, (C.17) 


(C.ll) 


+ 1 R + + 2 L 


\J 2u a/2 u a/2 
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FIG. 12. In order to realize a beam splitter with a non-chiral 
setup, two LLs with a local charging interaction which does 
not discriminate between the two wires is needed (indicicated 
by green oval, strength U > 0). In this setup, both intra-wire 
backscattering (short blue lines) and inter-wire backscattering 
(dashed red lines) are relevant, whereas inter-wire forward 
scattering (not shown in the figure) is irrelevant. The part 
of the wires inside the dashed boxed constitutes the quantum 
dot with a nonlocal charging interaction. 


and find for this tunneling operator the scaling dimension 

1 ( 1 ** 1 1 l\ 

01R.2L = j [Vi ^ + - + -j 

- \{Tu + 0 ’ <ai8) 

In the presence of a repulsive interaction U > 0 the 
charge mode velocity u > 1 , and hence giR^L < 1 makes 
interwire backscattering relevant. Similarly, we find 
that intra-wire backscattering is relevant with giR^iL = 
g\R,2L- Thus, this model allows to realize the intermedi¬ 
ate fixed point discussed in the main text in the frame¬ 
work of non-clirial LL wires. 


Appendix D: Solution of the refermionized model 
for g = 1 


In this appendix, we discuss the derivation of the scat¬ 
tering probabilities Eqs. (47a), (47b). 


1. Solution for a single scattering amplitude 


We first consider the simple situation where j v = 0 and 
only 7 h 7 ^ 0. In order to construct scattering states as 
eigenstates of the Hamiltonian, we need to make assump¬ 
tions about the commutators of 7747/4 with j\ and 4/i. If 
T/i 774 commutes with both fi and Ti, we will arrive at 
the standard solution of refermionization 32 . As there is 
no reason to expect that our result should differ from the 
standard one, we will make this assumption. Then, we 
define the new Majorana operator fi = * 7 / 17 / 4 / 1 , and with 


H = Hq + H tun find the following equations of motion 


id t ^i{x) = [H, 'I'i(x)] 


(D.la) 

= id x ^i{x) + a/ 27ra j h /1 5{x) 

id t ^\{x) = [.ff)'I'l(x)] 


(D.lb) 

= id x ^\(x) 

- a/2t ra 7 ^ f x S(x) 

-idth = [H, fi] 


(D.lc) 

= 2\Z‘2na 7 h 

*!(0)-*t(0) 



Away from x = 0, 4'i(x) satisfies the free fermion equa¬ 
tion and has plane waves as a solution. At x = 0, the field 
is discontinuous and acquires a phase shift. As a conse¬ 
quence, in the above equations of motion, 4/(0) needs to 
be interpreted as |['ki(0-|-)-|-’I'i(0—)]. We solve the equa¬ 
tions of motions by using the most general expression for 
scattering eigenstates 


4^1, e — J ^i + an j£ /i^ (D.2) 

For the wave functions, we make the ansatz 

<pn, e (x) = 9{-x)e i£X + 9(x)e iSll e i£X (D.3a) 
¥>i 4 i 6 (a:) = e{x)e iSli e i£X . (D.3b) 


Following Ref. 1321 we interpret the scattering states in the 
following way: the region x < 0 corresponds to the in¬ 
coming wave packet, the region x > 0 the scattered com¬ 
ponents of the wave packet. The wave function p>\\ mul¬ 
tiplying the creation operator ip* describes the amplitude 
for forward scattering, i.e. is associated with the compo¬ 
nent of the scattering state which continues to propagate 
along segment 1 of the edge, in the interior of the QD. 
The wave function multiplying the annihilation oper¬ 
ator %l>i on the other hand is interpreted as the amplitude 
for scattering onto the outgoing part of edge segment 4. 
Clearly, unitarity demands that the sum of the probabili¬ 
ties for forward and backwards scattering is equal to one, 
I Wn 1 2 + I Wi41 2 = 1- 


Using the ansatz Eq. (D.3b) in the equations of motion 
Eqs. ( D.la[ )-( |D.lc ), we find 


* [l — e tSl1 ] — 2v / 27ra7/ 1 aii = 0 
— ie lSli + 2v / 27ra7fta 1 i = 0 

eon = ^7 /,[- (l + e w “)+e*“] 

and finally the solutions 


ie — 47 ra 7 2 


(D.4a) 

(D.4b) 

(D.4c) 


(D.5a) 


e iSi4 = ~ 47ra 7 1 

ie — 47 ra 7 2 


(D.5b) 
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We note that, different from the case of truely non¬ 
interacting fermions, ||^ and |(/?i 4 j£ | 2 cannot be inter¬ 
preted as scattering probabilities for each value of energy 
separately. Instead, in order to obtain the scale depen¬ 
dent scattering probability, they need to be integrated 
over the derivative of the Fermi function to obtain 

Pl 4 = J dw(-/'(w))|(^i 4, w | 2 

2T h T , 

« ' arctan — . (D. 6 ) 

J- h 

In the last step, the derivative of the Fermi function was 
approximated as a box function of width temperature T 
and height 1/P. Identifying the width of the resonance as 
Tft, = 47 ra 7 ^, and introducing the scaling function F(x) = 
^ arctan |, this result can be written in a more compact 
form as Pi 4 = F(T/T). The scaling function has the 
limiting behaviors F(x) = 1 — \(x/2) 2 for small x and 
F(x ) = | for large x. Using these asymptotics, we find 

r ^ for T > T h 
P14 = < 2 (D. 7 ) 

[ 1 - for T < T h 

This result agrees with the standard perturbative analy¬ 
sis of a single scatterer in a LL. 

We now discuss a complementary scenario, in which 
7 /j = 0 and j v 7 ^ 0. By symmetry, it is clear that the 
solution can be found in an analogous way as discussed 
above. For future reference, the solution for the scatter¬ 
ing eigenstate is given by 

k,e = f (^22,6^2 + ^44,e^2 + 022,e/l) ■ (D.8) 


Now, for x < 0 the operator t/ 4 describes an incoming 
state along countour 2. For x > 0, the operator ip\ cre¬ 
ates a state which is forward scattered into the interior of 
the QD along contour 2, whereas i/j 2 describes a partial 
wave which is scattered onto the outgoing part of edge 
4. In what follows, it will be important that both tpi(x) 
and ip 2 (x) create an outgoing wave on edge 4. 

2. Solution for two competing scattering 
amplitudes 

We now consider the interesting case in which both 7 ^ 
and 7 t, are nonzero. Due to the Klein factors represented 
by products of Majorana fermions 774774 and 172174 , the 
two different scattering operators do not commute with 
each other, and the respective equations of motion are 
coupled. In addition to f± = 7771774 / 1 , we now introduce 
/2 = 7772774/1 ■ Now we need to take into account the 
nontrivial anti-commutator 

{/ 1 J 2 } = -277i/i7? 2 / 2 = -2C , (D.9) 

where we have made the assumption that /1 and f 2 com¬ 
mute with each other and with 771174 and 772774 . Since C 
is a product of an even number of Majoranas, it trivially 
commutes with T 1 , Tj, T 2 , and In addition, 

\fi,C\ = 0 , and [f 2 ,C] = 0 =► [H tun ,C]=0. 

(D.10) 

As C commutes with all operators, it commutes with the 
Hamiltonian and hence does not have any dynamics. In 
addition, it squares to one C 2 = 1. With this in mind, we 
find the following equations of motion for the full system: 


- id A hi Or) 
-idt^\(x) 
-idtfi 
-i d t '& 2 (x) 
-idt'f’lix) 

-idth 


[H, ^i(a:)] = idx'&xix) + V2ira f 1 8(x) 


[H, ^i(x)] = id x ^\(x) - V2ira 7 ^ /1 S(x) 


[H, / 1 ] = 2 V / 27ra [ 7 h*i(0) - 7fc*i(0) - 6 7 **2(0) + C lv ^>\( 0) 


[H, T 2 (a;)] = id x ^ 2 (x) + V2ira f 2 5(x) 


[H,%(x)] = id x \it(x) - V2na f 2 6(x) 


[H,f 2 ] = 2^a [ 7 ^ 2 ( 0 )- 7 ^ 2 ( 0 ) - C 7 ^i(0) + C 7^ t (0) 


(D.lla) 

(D.llb) 

(D.llc) 

(D.lld) 

(D.lle) 

(D.llf) 


In the following, we will derive scattering states as a so¬ 
lution of these equations by using the fact that C 2 = 1. 


We make the most general ansatz for operators 

= / (vil,«^I + ^ 14,€^1 + ¥ , 12 , e ^ , 2 + ^14, £^2 


+ail,e/l + «12,e/2 


(D.12a) 



T 3 22,e v k 2 + ^24 e *2 + < 7 3 21,e V f'i + P 2 4 e^l 
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creating eigenstates of the Hamiltonian H and satisfying 

eAl £ = [H,A{ e \, c4 e =[ff,4,e] ■ (D-13) 


Eqs. (D.15e D.15f), we can solve for 

ie — 47ra7^ 


e iSl1 = 


ie - 4?to (72 + 72 ) 


(D.16a) 


In the following, we only discuss the expression for A\ f . 

as the corresponding expression for A\ e is obtained by in¬ 
terchanging 7 /, with 7 „. Specifically, we make the ansatz 


ip\i,e(x) = 9(-x)e iex + 9(x)e iSll e iex 

(D.14a) 

Vu,e(x) = 0{x)e i5li e Ux 

(D.14b) 

Vi 2 ,t{x) = 9(x)e iSl2 e iex 

(D.14c) 

= 0{x)e^ x 

(D.14d) 


Imposing the condition Eq. (D.13), the equations 


i [1 - e iSl1 ] - 2 v / 2 ™ 


lh 


Oll,e — Cai2 e 


= 0 


— ie lSli + 2\ / '2na'yh ffln, e — Ca\ 2 , e =0 


— ie lSl2 — 2v / 27ra7„ 


«i2,e ^ Can,1 


= 0 


• £<5i 

-ze 1 


+ 2\/2iraj v ai 2 ,e — Can )e =0 


ea 11 ,e=^[-(l + e^)7. + e i5 -7/ l ] 


eai2, e — 


—e ,,5l2 7 w + e lSlir i v 


(D.15a) 

(D.15b) 

(D.15c) 

(D.15d) 

(D.15e) 

(D.15f) 


are found. Using C 2 = 1, we see that the first four 
equations only depend on the combination ttn, e — Cai 2 , e , 
and if we take an appropriate linear combination of 


e z5 14 = 

ie - 47ra(7^ + 7 2 ) 


(D.16b) 


e iSi 2 = -4tt ajh'Yy 

ie - 47ra(7 2 + 7 2 ) 


(D.16c) 




kKci'yh'Yv 


ie - 47ra(7 2 + 7 2 ) 


(D.16d) 


The probability for forward scattering is given by \e tSl1 1 2 , 
the probability for scattering onto lead 2 is given by 
|e i<5l2 | 2 , and the probability for scattering onto lead 4 
is given by |e *' Sl4 | 2 + \ e lSli \ 2 . There are two contributions 
for scattering onto lead 4, one due to the operator ipi in 
the scattering state Eq. (D.12a), and a second one due 


to the operator U 2 in the same scattering state. We will 
see in a moment that indeed both these contributions are 
needed in order to reproduce the perturbative result for 
the tunneling probability from lead 1 onto lead 4. 

Integrating over the derivative of the Fermi function 
as in the case of a single scatterer, we obtain the proba¬ 
bilities for transmission 



(D.17a) 

(D.17b) 

(D.17c) 


We now argue that these results agree with a perturba¬ 
tive calculation in 7 ^ and 7 v . In the limit of large T/T , 
one finds P 14 = 2 r^/T, exactly the same result as in 
the case j v = 0 discussed previously. This is to be exp- 
tected, since corrections due to 7 ,, can enter only in an 
additive fashion, and have to be of order 0 ( 7 ^ 7 ^). On the 
other hand, in the same limit, one finds P 12 = 2ThT v /TT, 
which is again to be expected since scattering from lead 
1 into lead 2 has to be a two-step process and has to be 
proportional to 7 ^ 7 ^. 
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